3.2.16 \(\int \frac {x}{\sqrt {a+b x+c x^2} (d+e x+f x^2)} \, dx\)

Optimal. Leaf size=402 \[ \frac {\left (e-\sqrt {e^2-4 d f}\right ) \tanh ^{-1}\left (\frac {4 a f+2 x \left (b f-c \left (e-\sqrt {e^2-4 d f}\right )\right )-b \left (e-\sqrt {e^2-4 d f}\right )}{2 \sqrt {2} \sqrt {a+b x+c x^2} \sqrt {2 a f^2-\sqrt {e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}\right )}{\sqrt {2} \sqrt {e^2-4 d f} \sqrt {2 a f^2-\sqrt {e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}-\frac {\left (\sqrt {e^2-4 d f}+e\right ) \tanh ^{-1}\left (\frac {4 a f+2 x \left (b f-c \left (\sqrt {e^2-4 d f}+e\right )\right )-b \left (\sqrt {e^2-4 d f}+e\right )}{2 \sqrt {2} \sqrt {a+b x+c x^2} \sqrt {2 a f^2+\sqrt {e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}\right )}{\sqrt {2} \sqrt {e^2-4 d f} \sqrt {2 a f^2+\sqrt {e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}} \]

________________________________________________________________________________________

Rubi [A]  time = 0.96, antiderivative size = 402, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 3, integrand size = 28, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.107, Rules used = {1032, 724, 206} \begin {gather*} \frac {\left (e-\sqrt {e^2-4 d f}\right ) \tanh ^{-1}\left (\frac {4 a f+2 x \left (b f-c \left (e-\sqrt {e^2-4 d f}\right )\right )-b \left (e-\sqrt {e^2-4 d f}\right )}{2 \sqrt {2} \sqrt {a+b x+c x^2} \sqrt {2 a f^2-\sqrt {e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}\right )}{\sqrt {2} \sqrt {e^2-4 d f} \sqrt {2 a f^2-\sqrt {e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}-\frac {\left (\sqrt {e^2-4 d f}+e\right ) \tanh ^{-1}\left (\frac {4 a f+2 x \left (b f-c \left (\sqrt {e^2-4 d f}+e\right )\right )-b \left (\sqrt {e^2-4 d f}+e\right )}{2 \sqrt {2} \sqrt {a+b x+c x^2} \sqrt {2 a f^2+\sqrt {e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}}\right )}{\sqrt {2} \sqrt {e^2-4 d f} \sqrt {2 a f^2+\sqrt {e^2-4 d f} (c e-b f)-b e f-2 c d f+c e^2}} \end {gather*}

Antiderivative was successfully verified.

[In]

Int[x/(Sqrt[a + b*x + c*x^2]*(d + e*x + f*x^2)),x]

[Out]

((e - Sqrt[e^2 - 4*d*f])*ArcTanh[(4*a*f - b*(e - Sqrt[e^2 - 4*d*f]) + 2*(b*f - c*(e - Sqrt[e^2 - 4*d*f]))*x)/(
2*Sqrt[2]*Sqrt[c*e^2 - 2*c*d*f - b*e*f + 2*a*f^2 - (c*e - b*f)*Sqrt[e^2 - 4*d*f]]*Sqrt[a + b*x + c*x^2])])/(Sq
rt[2]*Sqrt[e^2 - 4*d*f]*Sqrt[c*e^2 - 2*c*d*f - b*e*f + 2*a*f^2 - (c*e - b*f)*Sqrt[e^2 - 4*d*f]]) - ((e + Sqrt[
e^2 - 4*d*f])*ArcTanh[(4*a*f - b*(e + Sqrt[e^2 - 4*d*f]) + 2*(b*f - c*(e + Sqrt[e^2 - 4*d*f]))*x)/(2*Sqrt[2]*S
qrt[c*e^2 - 2*c*d*f - b*e*f + 2*a*f^2 + (c*e - b*f)*Sqrt[e^2 - 4*d*f]]*Sqrt[a + b*x + c*x^2])])/(Sqrt[2]*Sqrt[
e^2 - 4*d*f]*Sqrt[c*e^2 - 2*c*d*f - b*e*f + 2*a*f^2 + (c*e - b*f)*Sqrt[e^2 - 4*d*f]])

Rule 206

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1*ArcTanh[(Rt[-b, 2]*x)/Rt[a, 2]])/(Rt[a, 2]*Rt[-b, 2]), x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 724

Int[1/(((d_.) + (e_.)*(x_))*Sqrt[(a_.) + (b_.)*(x_) + (c_.)*(x_)^2]), x_Symbol] :> Dist[-2, Subst[Int[1/(4*c*d
^2 - 4*b*d*e + 4*a*e^2 - x^2), x], x, (2*a*e - b*d - (2*c*d - b*e)*x)/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a,
b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[2*c*d - b*e, 0]

Rule 1032

Int[((g_.) + (h_.)*(x_))/(((a_) + (b_.)*(x_) + (c_.)*(x_)^2)*Sqrt[(d_.) + (e_.)*(x_) + (f_.)*(x_)^2]), x_Symbo
l] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Dist[(2*c*g - h*(b - q))/q, Int[1/((b - q + 2*c*x)*Sqrt[d + e*x + f*x^2])
, x], x] - Dist[(2*c*g - h*(b + q))/q, Int[1/((b + q + 2*c*x)*Sqrt[d + e*x + f*x^2]), x], x]] /; FreeQ[{a, b,
c, d, e, f, g, h}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[e^2 - 4*d*f, 0] && PosQ[b^2 - 4*a*c]

Rubi steps

\begin {align*} \int \frac {x}{\sqrt {a+b x+c x^2} \left (d+e x+f x^2\right )} \, dx &=-\left (\left (-1-\frac {e}{\sqrt {e^2-4 d f}}\right ) \int \frac {1}{\left (e+\sqrt {e^2-4 d f}+2 f x\right ) \sqrt {a+b x+c x^2}} \, dx\right )+\left (1-\frac {e}{\sqrt {e^2-4 d f}}\right ) \int \frac {1}{\left (e-\sqrt {e^2-4 d f}+2 f x\right ) \sqrt {a+b x+c x^2}} \, dx\\ &=-\left (\left (2 \left (1-\frac {e}{\sqrt {e^2-4 d f}}\right )\right ) \operatorname {Subst}\left (\int \frac {1}{16 a f^2-8 b f \left (e-\sqrt {e^2-4 d f}\right )+4 c \left (e-\sqrt {e^2-4 d f}\right )^2-x^2} \, dx,x,\frac {4 a f-b \left (e-\sqrt {e^2-4 d f}\right )-\left (-2 b f+2 c \left (e-\sqrt {e^2-4 d f}\right )\right ) x}{\sqrt {a+b x+c x^2}}\right )\right )-\left (2 \left (1+\frac {e}{\sqrt {e^2-4 d f}}\right )\right ) \operatorname {Subst}\left (\int \frac {1}{16 a f^2-8 b f \left (e+\sqrt {e^2-4 d f}\right )+4 c \left (e+\sqrt {e^2-4 d f}\right )^2-x^2} \, dx,x,\frac {4 a f-b \left (e+\sqrt {e^2-4 d f}\right )-\left (-2 b f+2 c \left (e+\sqrt {e^2-4 d f}\right )\right ) x}{\sqrt {a+b x+c x^2}}\right )\\ &=-\frac {\left (1-\frac {e}{\sqrt {e^2-4 d f}}\right ) \tanh ^{-1}\left (\frac {4 a f-b \left (e-\sqrt {e^2-4 d f}\right )+2 \left (b f-c \left (e-\sqrt {e^2-4 d f}\right )\right ) x}{2 \sqrt {2} \sqrt {c e^2-2 c d f-b e f+2 a f^2-(c e-b f) \sqrt {e^2-4 d f}} \sqrt {a+b x+c x^2}}\right )}{\sqrt {2} \sqrt {c e^2-2 c d f-b e f+2 a f^2-(c e-b f) \sqrt {e^2-4 d f}}}-\frac {\left (1+\frac {e}{\sqrt {e^2-4 d f}}\right ) \tanh ^{-1}\left (\frac {4 a f-b \left (e+\sqrt {e^2-4 d f}\right )+2 \left (b f-c \left (e+\sqrt {e^2-4 d f}\right )\right ) x}{2 \sqrt {2} \sqrt {c e^2-2 c d f-b e f+2 a f^2+(c e-b f) \sqrt {e^2-4 d f}} \sqrt {a+b x+c x^2}}\right )}{\sqrt {2} \sqrt {c e^2-2 c d f-b e f+2 a f^2+(c e-b f) \sqrt {e^2-4 d f}}}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]  time = 0.97, size = 407, normalized size = 1.01 \begin {gather*} \frac {-\frac {\left (\sqrt {e^2-4 d f}+e\right ) \tanh ^{-1}\left (\frac {4 a f-b \left (\sqrt {e^2-4 d f}+e-2 f x\right )-2 c x \left (\sqrt {e^2-4 d f}+e\right )}{2 \sqrt {2} \sqrt {a+x (b+c x)} \sqrt {f \left (2 a f-b \left (\sqrt {e^2-4 d f}+e\right )\right )+c \left (e \sqrt {e^2-4 d f}-2 d f+e^2\right )}}\right )}{\sqrt {e^2-4 d f} \sqrt {f \left (2 a f-b \left (\sqrt {e^2-4 d f}+e\right )\right )+c \left (e \sqrt {e^2-4 d f}-2 d f+e^2\right )}}-\frac {\left (1-\frac {e}{\sqrt {e^2-4 d f}}\right ) \tanh ^{-1}\left (\frac {4 a f+b \left (\sqrt {e^2-4 d f}-e+2 f x\right )+2 c x \left (\sqrt {e^2-4 d f}-e\right )}{2 \sqrt {2} \sqrt {a+x (b+c x)} \sqrt {f \left (2 a f+b \left (\sqrt {e^2-4 d f}-e\right )\right )+c \left (-e \sqrt {e^2-4 d f}-2 d f+e^2\right )}}\right )}{\sqrt {f \left (2 a f+b \left (\sqrt {e^2-4 d f}-e\right )\right )+c \left (-e \sqrt {e^2-4 d f}-2 d f+e^2\right )}}}{\sqrt {2}} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[x/(Sqrt[a + b*x + c*x^2]*(d + e*x + f*x^2)),x]

[Out]

(-(((e + Sqrt[e^2 - 4*d*f])*ArcTanh[(4*a*f - 2*c*(e + Sqrt[e^2 - 4*d*f])*x - b*(e + Sqrt[e^2 - 4*d*f] - 2*f*x)
)/(2*Sqrt[2]*Sqrt[c*(e^2 - 2*d*f + e*Sqrt[e^2 - 4*d*f]) + f*(2*a*f - b*(e + Sqrt[e^2 - 4*d*f]))]*Sqrt[a + x*(b
 + c*x)])])/(Sqrt[e^2 - 4*d*f]*Sqrt[c*(e^2 - 2*d*f + e*Sqrt[e^2 - 4*d*f]) + f*(2*a*f - b*(e + Sqrt[e^2 - 4*d*f
]))])) - ((1 - e/Sqrt[e^2 - 4*d*f])*ArcTanh[(4*a*f + 2*c*(-e + Sqrt[e^2 - 4*d*f])*x + b*(-e + Sqrt[e^2 - 4*d*f
] + 2*f*x))/(2*Sqrt[2]*Sqrt[c*(e^2 - 2*d*f - e*Sqrt[e^2 - 4*d*f]) + f*(2*a*f + b*(-e + Sqrt[e^2 - 4*d*f]))]*Sq
rt[a + x*(b + c*x)])])/Sqrt[c*(e^2 - 2*d*f - e*Sqrt[e^2 - 4*d*f]) + f*(2*a*f + b*(-e + Sqrt[e^2 - 4*d*f]))])/S
qrt[2]

________________________________________________________________________________________

IntegrateAlgebraic [C]  time = 0.43, size = 204, normalized size = 0.51 \begin {gather*} \text {RootSum}\left [\text {$\#$1}^4 f-2 \text {$\#$1}^3 \sqrt {c} e-2 \text {$\#$1}^2 a f+\text {$\#$1}^2 b e+4 \text {$\#$1}^2 c d+2 \text {$\#$1} a \sqrt {c} e-4 \text {$\#$1} b \sqrt {c} d+a^2 f-a b e+b^2 d\&,\frac {\text {$\#$1}^2 \log \left (-\text {$\#$1}+\sqrt {a+b x+c x^2}-\sqrt {c} x\right )-a \log \left (-\text {$\#$1}+\sqrt {a+b x+c x^2}-\sqrt {c} x\right )}{2 \text {$\#$1}^3 f-3 \text {$\#$1}^2 \sqrt {c} e-2 \text {$\#$1} a f+\text {$\#$1} b e+4 \text {$\#$1} c d+a \sqrt {c} e-2 b \sqrt {c} d}\&\right ] \end {gather*}

Antiderivative was successfully verified.

[In]

IntegrateAlgebraic[x/(Sqrt[a + b*x + c*x^2]*(d + e*x + f*x^2)),x]

[Out]

RootSum[b^2*d - a*b*e + a^2*f - 4*b*Sqrt[c]*d*#1 + 2*a*Sqrt[c]*e*#1 + 4*c*d*#1^2 + b*e*#1^2 - 2*a*f*#1^2 - 2*S
qrt[c]*e*#1^3 + f*#1^4 & , (-(a*Log[-(Sqrt[c]*x) + Sqrt[a + b*x + c*x^2] - #1]) + Log[-(Sqrt[c]*x) + Sqrt[a +
b*x + c*x^2] - #1]*#1^2)/(-2*b*Sqrt[c]*d + a*Sqrt[c]*e + 4*c*d*#1 + b*e*#1 - 2*a*f*#1 - 3*Sqrt[c]*e*#1^2 + 2*f
*#1^3) & ]

________________________________________________________________________________________

fricas [B]  time = 14.50, size = 11311, normalized size = 28.14

result too large to display

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x/(c*x^2+b*x+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="fricas")

[Out]

1/4*sqrt(2)*sqrt((2*c*d^2 - b*d*e + a*e^2 - 2*a*d*f + (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*
b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt
((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5
+ (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^
4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*
a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3
 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (
b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*
f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e
^2)*f))*log(-(2*b^2*d^3 - 4*a*b*d^2*e + 2*a^2*d*e^2 + sqrt(2)*(b^2*d^2*e^2 - 2*a*b*d*e^3 + a^2*e^4 - 4*(b^2*d^
3 - 2*a*b*d^2*e + a^2*d*e^2)*f - (2*c^3*d^4*e^2 - 3*b*c^2*d^3*e^3 - 2*a*b*c*d*e^5 + a^2*c*e^6 + 8*a^3*d^2*f^4
+ (b^2*c + 3*a*c^2)*d^2*e^4 - 2*(2*a^2*b*d^2*e + 3*a^3*d*e^2 - 4*(a*b^2 - 3*a^2*c)*d^3)*f^3 + (5*a^2*b*d*e^3 +
 a^3*e^4 - 8*(b^2*c - 3*a*c^2)*d^4 + 4*(b^3 - 2*a*b*c)*d^3*e - 2*(5*a*b^2 - 11*a^2*c)*d^2*e^2)*f^2 - (8*c^3*d^
5 - 12*b*c^2*d^4*e + a^2*b*e^5 + 2*(b^2*c + 9*a*c^2)*d^3*e^2 + (b^3 - 10*a*b*c)*d^2*e^3 - 2*(a*b^2 - 4*a^2*c)*
d*e^4)*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6
- 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3
*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8
*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*
a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^
3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))*sqrt(c*x^2 + b*x + a)*sqrt((2*c
*d^2 - b*d*e + a*e^2 - 2*a*d*f + (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(
b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((b^2*d^2 - 2*a*b*d*e
 + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)
*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2
*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c
 - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c
)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^
2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^
2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)) + (4*b*c*d^3
+ a*b*d*e^2 - (b^2 + 4*a*c)*d^2*e)*x - (2*a*c^2*d^3*e^2 - 2*a*b*c*d^2*e^3 + 2*a^2*c*d*e^4 - 8*a^3*d^2*f^3 + 2*
(4*a^2*b*d^2*e + a^3*d*e^2 - 4*(a*b^2 - 2*a^2*c)*d^3)*f^2 - 2*(4*a*c^2*d^4 - 4*a*b*c*d^3*e + a^2*b*d*e^3 - (a*
b^2 - 6*a^2*c)*d^2*e^2)*f + (b*c^2*d^3*e^2 - b^2*c*d^2*e^3 + a*b*c*d*e^4 - 4*a^2*b*d^2*f^3 + (4*a*b^2*d^2*e +
a^2*b*d*e^2 - 4*(b^3 - 2*a*b*c)*d^3)*f^2 - (4*b*c^2*d^4 - 4*b^2*c*d^3*e + a*b^2*d*e^3 - (b^3 - 6*a*b*c)*d^2*e^
2)*f)*x)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 -
 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*
b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*
(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a
^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3
)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/x) - 1/4*sqrt(2)*sqrt((2*c*d^2 -
 b*d*e + a*e^2 - 2*a*d*f + (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 -
2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2
*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e
^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*
d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b
*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)
*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3
- 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2
- 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f))*log(-(2*b^2*d^3 - 4
*a*b*d^2*e + 2*a^2*d*e^2 - sqrt(2)*(b^2*d^2*e^2 - 2*a*b*d*e^3 + a^2*e^4 - 4*(b^2*d^3 - 2*a*b*d^2*e + a^2*d*e^2
)*f - (2*c^3*d^4*e^2 - 3*b*c^2*d^3*e^3 - 2*a*b*c*d*e^5 + a^2*c*e^6 + 8*a^3*d^2*f^4 + (b^2*c + 3*a*c^2)*d^2*e^4
 - 2*(2*a^2*b*d^2*e + 3*a^3*d*e^2 - 4*(a*b^2 - 3*a^2*c)*d^3)*f^3 + (5*a^2*b*d*e^3 + a^3*e^4 - 8*(b^2*c - 3*a*c
^2)*d^4 + 4*(b^3 - 2*a*b*c)*d^3*e - 2*(5*a*b^2 - 11*a^2*c)*d^2*e^2)*f^2 - (8*c^3*d^5 - 12*b*c^2*d^4*e + a^2*b*
e^5 + 2*(b^2*c + 9*a*c^2)*d^3*e^2 + (b^3 - 10*a*b*c)*d^2*e^3 - 2*(a*b^2 - 4*a^2*c)*d*e^4)*f)*sqrt((b^2*d^2 - 2
*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 +
2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c
 + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 -
8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 +
 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b
*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))*sqrt(c*x^2 + b*x + a)*sqrt((2*c*d^2 - b*d*e + a*e^2 - 2*a*
d*f + (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*
c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 -
2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a
^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*
b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 -
20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 -
 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*
c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2
)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)) + (4*b*c*d^3 + a*b*d*e^2 - (b^2 + 4*a*c)
*d^2*e)*x - (2*a*c^2*d^3*e^2 - 2*a*b*c*d^2*e^3 + 2*a^2*c*d*e^4 - 8*a^3*d^2*f^3 + 2*(4*a^2*b*d^2*e + a^3*d*e^2
- 4*(a*b^2 - 2*a^2*c)*d^3)*f^2 - 2*(4*a*c^2*d^4 - 4*a*b*c*d^3*e + a^2*b*d*e^3 - (a*b^2 - 6*a^2*c)*d^2*e^2)*f +
 (b*c^2*d^3*e^2 - b^2*c*d^2*e^3 + a*b*c*d*e^4 - 4*a^2*b*d^2*f^3 + (4*a*b^2*d^2*e + a^2*b*d*e^2 - 4*(b^3 - 2*a*
b*c)*d^3)*f^2 - (4*b*c^2*d^4 - 4*b^2*c*d^3*e + a*b^2*d*e^3 - (b^3 - 6*a*b*c)*d^2*e^2)*f)*x)*sqrt((b^2*d^2 - 2*
a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2
*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c
+ 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8
*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 +
2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*
c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/x) + 1/4*sqrt(2)*sqrt((2*c*d^2 - b*d*e + a*e^2 - 2*a*d*f -
(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^
3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^
3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2
 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d
^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b
^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c
^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d
*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2
- (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f))*log(-(2*b^2*d^3 - 4*a*b*d^2*e + 2*a^2*d*e^2 +
sqrt(2)*(b^2*d^2*e^2 - 2*a*b*d*e^3 + a^2*e^4 - 4*(b^2*d^3 - 2*a*b*d^2*e + a^2*d*e^2)*f + (2*c^3*d^4*e^2 - 3*b*
c^2*d^3*e^3 - 2*a*b*c*d*e^5 + a^2*c*e^6 + 8*a^3*d^2*f^4 + (b^2*c + 3*a*c^2)*d^2*e^4 - 2*(2*a^2*b*d^2*e + 3*a^3
*d*e^2 - 4*(a*b^2 - 3*a^2*c)*d^3)*f^3 + (5*a^2*b*d*e^3 + a^3*e^4 - 8*(b^2*c - 3*a*c^2)*d^4 + 4*(b^3 - 2*a*b*c)
*d^3*e - 2*(5*a*b^2 - 11*a^2*c)*d^2*e^2)*f^2 - (8*c^3*d^5 - 12*b*c^2*d^4*e + a^2*b*e^5 + 2*(b^2*c + 9*a*c^2)*d
^3*e^2 + (b^3 - 10*a*b*c)*d^2*e^3 - 2*(a*b^2 - 4*a^2*c)*d*e^4)*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^
4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b
*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^
3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e -
 (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c
^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c
- 2*a^2*c^2)*d*e^4)*f)))*sqrt(c*x^2 + b*x + a)*sqrt((2*c*d^2 - b*d*e + a*e^2 - 2*a*d*f - (c^2*d^2*e^2 - b*c*d*
e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b
*e^3 - (b^2 - 6*a*c)*d*e^2)*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2
*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3
*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a
^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^
2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^
5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e
^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*
d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f)) + (4*b*c*d^3 + a*b*d*e^2 - (b^2 + 4*a*c)*d^2*e)*x + (2*a*c^2*d^3*e^
2 - 2*a*b*c*d^2*e^3 + 2*a^2*c*d*e^4 - 8*a^3*d^2*f^3 + 2*(4*a^2*b*d^2*e + a^3*d*e^2 - 4*(a*b^2 - 2*a^2*c)*d^3)*
f^2 - 2*(4*a*c^2*d^4 - 4*a*b*c*d^3*e + a^2*b*d*e^3 - (a*b^2 - 6*a^2*c)*d^2*e^2)*f + (b*c^2*d^3*e^2 - b^2*c*d^2
*e^3 + a*b*c*d*e^4 - 4*a^2*b*d^2*f^3 + (4*a*b^2*d^2*e + a^2*b*d*e^2 - 4*(b^3 - 2*a*b*c)*d^3)*f^2 - (4*b*c^2*d^
4 - 4*b^2*c*d^3*e + a*b^2*d*e^3 - (b^3 - 6*a*b*c)*d^2*e^2)*f)*x)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4
*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*
d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3
 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e -
(b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^
4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c -
 2*a^2*c^2)*d*e^4)*f)))/x) - 1/4*sqrt(2)*sqrt((2*c*d^2 - b*d*e + a*e^2 - 2*a*d*f - (c^2*d^2*e^2 - b*c*d*e^3 +
a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 -
 (b^2 - 6*a*c)*d*e^2)*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5
 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^
2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*
d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2
+ 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b
^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b
*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e
+ a*b*e^3 - (b^2 - 6*a*c)*d*e^2)*f))*log(-(2*b^2*d^3 - 4*a*b*d^2*e + 2*a^2*d*e^2 - sqrt(2)*(b^2*d^2*e^2 - 2*a*
b*d*e^3 + a^2*e^4 - 4*(b^2*d^3 - 2*a*b*d^2*e + a^2*d*e^2)*f + (2*c^3*d^4*e^2 - 3*b*c^2*d^3*e^3 - 2*a*b*c*d*e^5
 + a^2*c*e^6 + 8*a^3*d^2*f^4 + (b^2*c + 3*a*c^2)*d^2*e^4 - 2*(2*a^2*b*d^2*e + 3*a^3*d*e^2 - 4*(a*b^2 - 3*a^2*c
)*d^3)*f^3 + (5*a^2*b*d*e^3 + a^3*e^4 - 8*(b^2*c - 3*a*c^2)*d^4 + 4*(b^3 - 2*a*b*c)*d^3*e - 2*(5*a*b^2 - 11*a^
2*c)*d^2*e^2)*f^2 - (8*c^3*d^5 - 12*b*c^2*d^4*e + a^2*b*e^5 + 2*(b^2*c + 9*a*c^2)*d^3*e^2 + (b^3 - 10*a*b*c)*d
^2*e^3 - 2*(a*b^2 - 4*a^2*c)*d*e^4)*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2
*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2
 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b
^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2
*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^
2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))*sq
rt(c*x^2 + b*x + a)*sqrt((2*c*d^2 - b*d*e + a*e^2 - 2*a*d*f - (c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 - 4*a^2*d*f^3
 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*a*c)*d*e^2)
*f)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^
4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2 - 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3
 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*
c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*
c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3
*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/(c^2*d^2*e^2 - b*c*d*e^3 + a*c*e^4 -
4*a^2*d*f^3 + (4*a*b*d*e + a^2*e^2 - 4*(b^2 - 2*a*c)*d^2)*f^2 - (4*c^2*d^3 - 4*b*c*d^2*e + a*b*e^3 - (b^2 - 6*
a*c)*d*e^2)*f)) + (4*b*c*d^3 + a*b*d*e^2 - (b^2 + 4*a*c)*d^2*e)*x + (2*a*c^2*d^3*e^2 - 2*a*b*c*d^2*e^3 + 2*a^2
*c*d*e^4 - 8*a^3*d^2*f^3 + 2*(4*a^2*b*d^2*e + a^3*d*e^2 - 4*(a*b^2 - 2*a^2*c)*d^3)*f^2 - 2*(4*a*c^2*d^4 - 4*a*
b*c*d^3*e + a^2*b*d*e^3 - (a*b^2 - 6*a^2*c)*d^2*e^2)*f + (b*c^2*d^3*e^2 - b^2*c*d^2*e^3 + a*b*c*d*e^4 - 4*a^2*
b*d^2*f^3 + (4*a*b^2*d^2*e + a^2*b*d*e^2 - 4*(b^3 - 2*a*b*c)*d^3)*f^2 - (4*b*c^2*d^4 - 4*b^2*c*d^3*e + a*b^2*d
*e^3 - (b^3 - 6*a*b*c)*d^2*e^2)*f)*x)*sqrt((b^2*d^2 - 2*a*b*d*e + a^2*e^2)/(c^4*d^4*e^2 - 2*b*c^3*d^3*e^3 - 2*
a*b*c^2*d*e^5 + a^2*c^2*e^6 - 4*a^4*d*f^5 + (b^2*c^2 + 2*a*c^3)*d^2*e^4 + (8*a^3*b*d*e + a^4*e^2 - 8*(a^2*b^2
- 2*a^3*c)*d^2)*f^4 - 2*(a^3*b*e^3 + 2*(b^4 - 4*a*b^2*c + 6*a^2*c^2)*d^3 - 4*(a*b^3 - a^2*b*c)*d^2*e + (a^2*b^
2 + 6*a^3*c)*d*e^2)*f^3 - (8*(b^2*c^2 - 2*a*c^3)*d^4 - 8*(b^3*c - a*b*c^2)*d^3*e - (b^4 - 20*a*b^2*c + 22*a^2*
c^2)*d^2*e^2 + 2*(a*b^3 - 5*a^2*b*c)*d*e^3 - (a^2*b^2 + 2*a^3*c)*e^4)*f^2 - 2*(2*c^4*d^5 - 4*b*c^3*d^4*e + a^2
*b*c*e^5 + (b^2*c^2 + 6*a*c^3)*d^3*e^2 + (b^3*c - 5*a*b*c^2)*d^2*e^3 - 2*(a*b^2*c - 2*a^2*c^2)*d*e^4)*f)))/x)

________________________________________________________________________________________

giac [F(-1)]  time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Timed out} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x/(c*x^2+b*x+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="giac")

[Out]

Timed out

________________________________________________________________________________________

maple [B]  time = 0.01, size = 1516, normalized size = 3.77

result too large to display

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x/(c*x^2+b*x+a)^(1/2)/(f*x^2+e*x+d),x)

[Out]

-1/2/(-4*d*f+e^2)^(1/2)/f*2^(1/2)/((2*a*f^2-b*e*f-2*c*d*f+c*e^2-(-4*d*f+e^2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*c*e)
/f^2)^(1/2)*ln(((b*f-c*e-(-4*d*f+e^2)^(1/2)*c)*(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f)/f+(2*a*f^2-b*e*f-2*c*d*f+c*e^2
-(-4*d*f+e^2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*c*e)/f^2+1/2*2^(1/2)*((2*a*f^2-b*e*f-2*c*d*f+c*e^2-(-4*d*f+e^2)^(1/
2)*b*f+(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/2)*(4*(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f)^2*c+4*(b*f-c*e-(-4*d*f+e^2)^(1/2
)*c)*(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f)/f+2*(2*a*f^2-b*e*f-2*c*d*f+c*e^2-(-4*d*f+e^2)^(1/2)*b*f+(-4*d*f+e^2)^(1/
2)*c*e)/f^2)^(1/2))/(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f))*e-1/2/f*2^(1/2)/((2*a*f^2-b*e*f-2*c*d*f+c*e^2-(-4*d*f+e^
2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/2)*ln(((b*f-c*e-(-4*d*f+e^2)^(1/2)*c)*(x+1/2*(e+(-4*d*f+e^2)^(1/2
))/f)/f+(2*a*f^2-b*e*f-2*c*d*f+c*e^2-(-4*d*f+e^2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*c*e)/f^2+1/2*2^(1/2)*((2*a*f^2-
b*e*f-2*c*d*f+c*e^2-(-4*d*f+e^2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/2)*(4*(x+1/2*(e+(-4*d*f+e^2)^(1/2))
/f)^2*c+4*(b*f-c*e-(-4*d*f+e^2)^(1/2)*c)*(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f)/f+2*(2*a*f^2-b*e*f-2*c*d*f+c*e^2-(-4
*d*f+e^2)^(1/2)*b*f+(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/2))/(x+1/2*(e+(-4*d*f+e^2)^(1/2))/f))-1/2/f*2^(1/2)/((2*a*
f^2-b*e*f-2*c*d*f+c*e^2+(-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/2)*ln(((b*f-c*e+(-4*d*f+e^2)^(1
/2)*c)*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)/f+(2*a*f^2-b*e*f-2*c*d*f+c*e^2+(-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1
/2)*c*e)/f^2+1/2*2^(1/2)*((2*a*f^2-b*e*f-2*c*d*f+c*e^2+(-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/
2)*(4*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)^2*c+4*(b*f-c*e+(-4*d*f+e^2)^(1/2)*c)*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)
/f+2*(2*a*f^2-b*e*f-2*c*d*f+c*e^2+(-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/2))/(x-1/2*(-e+(-4*d*
f+e^2)^(1/2))/f))+1/2/(-4*d*f+e^2)^(1/2)/f*2^(1/2)/((2*a*f^2-b*e*f-2*c*d*f+c*e^2+(-4*d*f+e^2)^(1/2)*b*f-(-4*d*
f+e^2)^(1/2)*c*e)/f^2)^(1/2)*ln(((b*f-c*e+(-4*d*f+e^2)^(1/2)*c)*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)/f+(2*a*f^2-b
*e*f-2*c*d*f+c*e^2+(-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e)/f^2+1/2*2^(1/2)*((2*a*f^2-b*e*f-2*c*d*f+c*e^
2+(-4*d*f+e^2)^(1/2)*b*f-(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/2)*(4*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)^2*c+4*(b*f-c*
e+(-4*d*f+e^2)^(1/2)*c)*(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f)/f+2*(2*a*f^2-b*e*f-2*c*d*f+c*e^2+(-4*d*f+e^2)^(1/2)*
b*f-(-4*d*f+e^2)^(1/2)*c*e)/f^2)^(1/2))/(x-1/2*(-e+(-4*d*f+e^2)^(1/2))/f))*e

________________________________________________________________________________________

maxima [F(-2)]  time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Exception raised: ValueError} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x/(c*x^2+b*x+a)^(1/2)/(f*x^2+e*x+d),x, algorithm="maxima")

[Out]

Exception raised: ValueError >> Computation failed since Maxima requested additional constraints; using the 'a
ssume' command before evaluation *may* help (example of legal syntax is 'assume(4*d*f-e^2>0)', see `assume?` f
or more details)Is 4*d*f-e^2 positive, negative or zero?

________________________________________________________________________________________

mupad [F]  time = 0.00, size = -1, normalized size = -0.00 \begin {gather*} \int \frac {x}{\sqrt {c\,x^2+b\,x+a}\,\left (f\,x^2+e\,x+d\right )} \,d x \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x/((a + b*x + c*x^2)^(1/2)*(d + e*x + f*x^2)),x)

[Out]

int(x/((a + b*x + c*x^2)^(1/2)*(d + e*x + f*x^2)), x)

________________________________________________________________________________________

sympy [F]  time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \int \frac {x}{\sqrt {a + b x + c x^{2}} \left (d + e x + f x^{2}\right )}\, dx \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x/(c*x**2+b*x+a)**(1/2)/(f*x**2+e*x+d),x)

[Out]

Integral(x/(sqrt(a + b*x + c*x**2)*(d + e*x + f*x**2)), x)

________________________________________________________________________________________